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We study three-terminal linear and nonlinear transport through an Aharonov-Bohm interferome- 
ter containing a double quantum dot using the nonequilibrium Green's function method. Under the 
condition that one of the three terminals is a voltage probe, we show that the linear conductance is 
symmetric with respect to the magnetic field (phase symmetry). However, in the nonlinear trans- 
port regime, the phase symmetry is broken. Unlike two-terminal transport, the phase symmetry is 
broken even in noninteracting electron systems. Based on the lowest-order nonlinear conductance 
coefficient with respect to the source-drain bias voltage, we discuss the direction in which the phase 
shifts with the magnetic field. When the higher harmonic components of the Aharonov-Bohm os- 
cillations are negligible, the phaseshift is a monotonically increasing function with respect to the 
source-drain bias voltage. To observe the Aharonov-Bohm oscillations with higher visibility, we need 
strong coupling between the quantum dots and the voltage probe. However, this leads to dephasing 
since the voltage probe acts as a Biittiker dephasing probe. The interplay between such antithetic 
concepts provides a peak in the visibility of the Aharonov-Bohm oscillations when the coupling 
between the quantum dots and the voltage probe changes. 

PACS numbers: 73.23.-b, 73.63.Kv, 73.40.Gk, 05.60.Gg 

I. INTRODUCTION 

The manifestation of quantum phase coherence forms one of the foundations of the physics of mesoscopic systems*, 
and is attracting the attention of many physicists. Quantum phase coherence is detectable by quantum interfer- 
ence experiments employing, for example, the Aharonov-Bohm (AB) effect^. In an AB interferometer containing a 
quantum dot (QD), the AB effects in the transport properties have been widely studied both theoretically^ and 
experimentally^—. The experiments show that phase coherence is maintained during the tunneling process through a 
QD. 

By applying a finite bias voltage across a QD system, we can easily realize a nonequilibrium steady state condition. 
Therefore, QD systems provide a suitable stage for testing theories related to nonequilibirum systems. It is reasonable 
to expect that driving a system out of equilibrium will provide a new understanding of quantum interference effects. 
Some of symmetries present at equilibrium, which underline a linear response, may be broken, and at the same time 
new qualitative features may emerge. 

The Onsager-Casimir symmetry relation states that the linear conductance in two-terminal systems should be 
symmetric with respect to an external magnetic field (phase rigidity )2ii£. In the nonlinear transport regime, however, 
it is not necessary for this phase symmetry to be satisfied. Recently, phase symmetry breaking in the nonlinear 
transport regime in two-terminal systems has been extensively studied both theoretically^^— and experimentally^— . 
Phase rigidity is not enforced in a two-terminal conductor if the conductor is interacting with another subsystem in 
a nonequilibrium situation 2 ^— . Moreover, in a multi-terminal conductor that includes the lossy channels, the phase 
symmetry breaks since the additional reservoirs allow losses of current and lead to the violation of unitarit y 24 ' 25 . 
In Ref. |26| . Biittiker had shown that the phase symmetry relation is satisfied in linear transport regime through 
multi-terminal device that all terminals except the source and drain reservoirs are voltage probes. In this paper, we 
consider a three-terminal system that satisfies unitarity where one of the three terminals is a voltage probe. The 
question is now, when the voltage probe is included in the AB interferometer, is the phase symmetry broken due to 
the voltage (electrochemical potential) fluctuation of the voltage probe in the nonlinear transport regime? The voltage 
probe is mathematically equivalent to the Biittiker dephasing probe 2 ^. In this approach, a system is connected to 
a virtual electron reservoir through a fictitious voltage probe. Electrons are scattered into such a probe, lose their 
phase memories with a certain probability, and are then reinjected into the system. Thus, the voltage probe induces 
dephasing and at the same time it assists quantum phase coherence as part of the AB interferometer. Moreover, in 
two-terminal systems, the phase symmetry is not broken in noninteracting electron systems. A voltage probe is an 
infinite impedance terminal with zero net current and imposes a constraint. Then, we expect that the phase symmetry 
may be broken due to this constraint even in non-interacting electron systems. 

In this paper, we study the phase and amplitude of the AB oscillation in a three-terminal AB interferometer 
device to address the following two main issues: (i) Is the phase symmetry broken in nonlinear transport through 
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a three-terminal AB interferometer that includes a voltage probe? To clarify this, we investigate the phase of the 
AB oscillation in the lowest-order nonlinear conductance coefficient with respect to the bias voltage, (ii) The other 
major issue this paper addresses is the dephasing effects caused by the voltage probe, which is a component of an AB 
interferometer. We examine the amplitude of AB oscillations in transport properties to discuss the way in which the 
interference effect is suppressed by coupling with the voltage probe. 

Here we investigate linear and nonlinear three-terminal transport through an AB interferometer containing a double 
quantum dot (DQD) using the nonequilibrium Green's function metho d 28 ' 29 . Recently, the AB effects in an AB 
interferometer containing a DQD have been thoroughly examined^r— . In laterally coupled DQD systems, coherent 
indirect coupling between two QDs via a reservoir is essential in terms of coherent transport through a DQ D 30 i 31 . 
We introduce the coherent indirect coupling parameter a, which characterizes the strength of the indirect coupling 
between two QDs via a reservoir. We consider three reservoirs, namely a source (S), a drain (D), and a voltage 
probe (tp) as shown in Fig. [TJ The electrochemical potential of the voltage probe is determined in order to satisfy 
the condition that the net current through the voltage probe vanishes. This is equivalent to a Biittiker dephasing 
probed. Thus, the coupling between the QDs and the voltage probe gives rise to the dephasing in the electronic states 
of the DQD. The coherent indirect coupling in a laterally coupled double quantum dot characterizes the coherence 
between two quantum dots. Thus, to study the interplay between the inter-dot coherence due to the coherent indirect 
coupling and the dephasing due to the voltage probe, the three-terminal system with a voltage probe coupled to 
a DQD is intersting. We show that the linear conductance is symmetric with respect to the magnetic field and 
show that this phase symmetry is broken in the nonlinear transport regime. In two-terminal systems, an electron- 
electron interaction is essential for breaking the phase symmetryiir— . However, in a three-terminal system including a 
voltage probe, we show that the phase symmetry is broken even in noninteracting electron systems. When the higher 
harmonic components of the AB oscillations are negligible, we derive an expression for the phaseshift and show that 
the phaseshift is independent of the coherent indirect coupling and monotonically increasing function with respect to 
the source-drain bias voltage under low source-drain bias voltage conditions. To observe AB oscillations with a large 
amplitude, we need strong coupling between the QD and the voltage probe. However, this induces the dephasing of 
the electronic states in the DQD. The competition between such antithetic concepts generates a peak structure when 
the coupling between the quantum dots and the voltage probe changes. 

The outline of this paper is as follows. In Sec. [ITJ we introduce a microscopic model Hamiltonian with the three 
reservoirs (S, D, and <p as shown in Fig. [TJ) and the notion of coherent indirect couplin g 30 ' 31 . In Sec. IIIIl we provide 
a theoretical formulation based on the nonequilibrium Green's function metho d 28 ' 29 . In particular, we impose a 
condition for the voltage probe (p. Section llVI is devoted to theoretical results for nonlinear transport properties. In 
Sec. El we discuss the interplay between the coherent effect of the coupling r v and the dephasing in relation to the 
Biittiker probe, where T v is the coupling strength between the QD and the voltage probe. Section [VTl summarizes our 
results. In Appendix^ we derive the relation that the transmission probability has to satisfy. We show that the linear 
conductance satisfies the phase symmetry relation in Appendix [B] In Appendix [C] we show the antisymmetricity of 
the lowest-order nonlinear conductance coefficient when the system has mirror symmetry. In Appendix [D] wc provide 
the detailed derivation of the visibility of the AB oscillations in the linear conductance in the limit of a = 1 and 

r„ -> oo. 



We consider an AB interferometer containing a DQD coupled to three reservoirs as shown in Fig. [TJ To focus on 
the coherent charge transport, we neglect the spin degree of freedom. Moreover, we assume that the level spacing is 
much larger than the source-drain bias voltage, and consider only a single energy level in each QD. The Hamiltonian 
represents the sum of the following terms: H = Hr + Hdqd + Ht- The Hamiltonian of the Fermi liquid reservoirs is 



where e u k is the electron energy with wave number k in the reservoir i/, and the operator a u k (a v k^) annihilates 
(creates) an electron in the reservoir v. Hdqd describes the isolated DQD, 



II. MODEL 
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FIG. 1: Schematic diagram of an AB interferometer containing a DQD coupled to three reservoirs. The two QDs couple to 
three reservoirs, namely the source (S), drain (D), and voltage probe (<p). S12 is the propagation length of electrons in the 
voltage probe. $ is a magnetic flux threading through the AB interferometer. The dashed line indicates the axis of the mirror 
symmetry. 



Here Ej is the energy level of the jth QD, and t c is the direct inter-dot tunnel coupling. The tunneling Hamiltonian 
between the QDs and the reservoirs is given by 

Ht = Y. [4fc°s* tc l + t[ Dk a Dk ] c 2 + tf k e^/ 2 a^ Cl + tf k e-^ 2 a^c 2 + h. 
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where is the tunneling amplitude between the jth QD and the reservoir v. As an effect of the magnetic flux, we 
introduced the Peierls phase factors e ±l */ 2 (</> = 27r$/$o is an AB phase, where <J> is the magnetic flux threading 
through an AB ring consisting of the two QDs and the voltage probe as shown in Fig. [TJ and $q = h/e is the magnetic 
flux quantum.) 

For the reservoir v G S, D, the line width functions are 
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which is assumed to be independent of the energy in the range of interest (p u is the density of states in the reservoir 
v), and when v G S(D), we have j = 1(2). Here we introduced the notation M^, which denotes the matrix 
element of a 2 x 2 matrix M , where the boldface notation indicates a 2 x 2 matrix whose basis is a localized state in 
each QD. In contrast, for the reservoir tp, the linewidth function matrix is not diagonal as follows 
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where we assumed a wide-band limit, namely we neglected the energy dependence. Using the matrix representation, 
we have 
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and the total lincwidth function matrix is defined as F(<f>) — F s + T D + T v (<fi). When we calculate the linewidth 

functions from the definitions (4) and (5), we estimate the tunneling amplitude t[5 in the tunneling Hamiltonian 
(3) using the Bardeen's formula 3 ^. In the Bardecn's theory, the tunneling amplitude can be expressed by the wave 
functions of evanescent mode of the reservoir v and a localized electron in the jth QD. Then, the tunneling amplitude 
^vl depends on the position of the quantum dot. As a result, the coherent indirect coupling parameter a is a function 
of the distance between two QDs. This coherent indirect coupling parameter a characterizes the strength of the 
indirect coupling between two QDs via the voltage prob o 30 i 31 . The coherent indirect coupling parameter becomes 
small and changes its sign with increasing distance (si2 in Fig. [T]) between the two QDs^i. The importance of the sign 
of the coherent indirect coupling parameter was pointed out by S. A. Gurvita 3 ^. The influence of the AB effects on 
the sign of the coherent indirect coupling parameters has been examined experimentally 3 ^. From Eq. all physical 
quantities are invariant under the transformation that we change the sign of a and shift the AB phase by tt. 



III. FORMULATION 



The tunneling current from the reservoir S to the DQD is given by^- 

Idqd{4>) = I J de[fs{e)-f D {e)]Tv{G r {e,<p)T s G a {e,<p)T D } 

J de[f s (e)-f v (e)]Tr{G r (e^)T s G a (e^)T*^)} 



e 
h 



= T / de[f s (e)-f D (e)]T DS (e, <!>) + - / de[f s (e) - f^T^i 



where the retarded Green's function is the Fourier transform of 

g^(m') = -ie(t-t')({ Ci (t), Cj Ht')}), 



(7) 



(8) 



and the advanced Green's function is obtained from the retarded Green's function: G a (e,<p) = [G r (e, (f>)]^ ■ /i/(e) is 
the Fermi-Dirac distribution function of the reservoir v defined as 
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where /i„ is the electrochemical potential of the reservoir v, and T is the temperature. In the following discussions, we 
assume that the source and drain reservoirs have electrochemical potentials /is = /i + eVsD/2 and fijj = fi— eVs£>/2 
with the source-drain bias voltage Vsd, and fi = 0. Here we define the transmission probability from the reservoir v 
to the reservoir £ as 



r c „(e, </>) = Tr {G r (e, 4>)T" '{<f)G a (e, <t>)T*{4>)} • 



(10) 



The electrochemical potential \i v is determined by the condition that the net current I v (4>) flowing through the voltage 
probe ip vanishes. Then, we impose following condition to determine 



dc[/ v (e, </>, V SD ) - fs(e)]T Sv (e, <£) + ^ / de[f v {e, 0, V SD ) - f D {e)]T Dlp {e, 



(11) 



The reservoir ip that satisfies such a condition is equivalent to the Biittiker dephasing reservoir—. 

To calculate the above physical quantities, we need the retarded Green's function. In our model, the retarded 
Green's function is given by 
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where 



A(e,4>) = 



e-ei 



~JT + 2 r22 



(12) 



(13) 



5 



In the following, we calculate the linear and the lowest-order nonlinear conductance coefficient with respect to the 
source-drain bias voltage. In general, the current is expressed as a polynomial function of the source-drain bias voltage 

V S D, 



IdqdW = G% D ^)V SD + -GWqdMVsb 
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To calculate the linear and nonlinear conductance coefficient, we focus on the zero temperature condition and employ 
the following linear approximation for the transmission probability 
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From the condition , the electrochemical potential of the reservoir ip is given by 
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Using Eqs. (J7)), (|14p. (fT5|) . and (|16[) . the linear conductance is given by 
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FIG. 2: AB oscillations of the linear conductance under the condition of mirror symmetry for V — T^, t c /KT = — 1. (a) 
td/KT = 1. (b) ed/hT — 0.2. The solid and broken lines indicate cases where a — 0.5 and a = —0.5, respectively. 



We confirmed the relation 

T iv {e,-4>)=T vi {e,4>). (25) 

This is based on time-reversal symmetry. In the discussions in the later sections, we sometimes consider a highly 
symmetric situation, namely ei = 62 = e c i, = Tjj = T, and T v i = T V 2 = ^u> (mirror symmetry with respect to the 
dashed line in Fig. [TJ. Under this condition, we can have following additional relation: 

T vS (e,<t>)=T Dv (e,4>). (26) 

Before we discuss nonlinear transport in the next section, we consider linear transport. From the conservation of 
probability, we have the following relation 

^T^(6,0) = ^T e „(e». (27) 

The derivation of this relation is given in Appendix |Al As proven in Appendix [Bj using this relation for arbitrary 
parameters, we find that the linear conductance satisfies the Onsager-Casimir symmetry relation^iiS 

G% D {-<i>)=G% D {<i>). (28) 

This satisfies the Buttiker's result for the linear conductance^. 

In Fig. we show the numerical results for the AB oscillations of the linear conductance under the condition of 
mirror symmetry. In Fig. H](a), we plot the AB oscillations of the linear conductance when T = T v , t c /HT = — 1, 
and €d/hT = 1. In this case, the convex shape of the linear conductance at (f> = depends on the sign of the coherent 
indirect coupling parameter a. Similarly, we show the AB oscillations of the linear conductance when T = T v , 
t c /hT = —1, and ed/Sr = 0.2. Under this condition, we find that the convex shape of the linear conductance at 4> = 
is independent of the sign of the coherent indirect coupling parameter a. 



IV. PHASE SYMMETRY BREAKING 



In this section, we discuss the phase symmetry breaking in the nonlinear transport regime. We derive the condition 
under which the phase symmetry is broken by calculating the lowest-order nonlinear conductance coefficient. 



A. Nonlinear transport in weak inter-dot coherent coupling 



Here we discuss the nonlinear transport under a finite source-drain bias voltage. Before discussing the general 
properties of the lowest-order nonlinear conductance coefficient, we focus on the weak inter-dot coupling situation 
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where |a|fiF v , \t c \ <C HT to remove the higher harmonic components of the AB oscillations and obtain an intuitive 
picture of the phaseshift. Under the mirror symmetry condition, the tunneling current through a DQD is written as 
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Here A<f> (— it < A<f> < w) satisfies the relations 
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Thus, we can directly derive the expression of the phaseshift from Eq. (I35[) 
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When Vsd ^ 0, the phaseshift is finite and the phase symmetry is broken. From Eq. (|37|) we find that the phaseshift 
is independent of the inter-dot couplings, namely t c and ahT^, and a monotonic function with respect to the source- 
drain bias voltage Vsd- In this expression, the phaseshift is defined in the range of — \ < A<fi < However, the 
original A<p is — n < A(f> < 7r range. When I\ is negative, namely ^ < A<p < tt or — tt < A<p < — -|, we can express 
the results in the — -| < A<f> < ^ range by shifting the phaseshift by n and changing the sign of the current. 





(a) 



(b) 



FIG. 3: QD energy e d dependence of the phaseshift for eV SD /hT = 0.1. (a) Phaseshift for = V (hT < 3hT v ). (b) 

v v = o.ir (hv > 3hT v ). 
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Here we consider the QD energy dependence of the phaseshift. First we consider the case where HT < 3fir v . We 
have 

lim Ac/> = ±^, 
and thus we find that the phaseshift jumps by ir at 

, ^{ZhT v -hT)(hT + hv v ) 

e d ,± = ± 

as shown in Fig. [3] (a) when = T and eVsu/^r = 0.1. For A = 0, the phaseshift is not defined since the linear 
conductance does not exhibit an Aharonov-Bohm (AB) oscillation as shown in Fig. Hta) dashed-line. Although the 
phaseshift jumps by 7r at A = 0, the tunneling current changes smoothly as shown in Fig. |4](b). This phaseshift jump 
is based on the fact that the convex shape of the linear conductance changes at A = (see Fig. |4] (a)). 

Next we consider the case where KT > ihT v . We find that the phaseshift is always positive, a smooth function 
with respect to the QD energy e^, and has its maximum value at = 0. Fig. [3](b) shows the numerical result when 
= O.ir and eVso/fiT = 0.1. 
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FIG. 4: AB oscillations of the linear conductance and the tunneling current through a DQD where A > (T^ = 0.99F), A = 
(1% = T), and A < (T v = l.Oir) when e d /hY = 1, t c /hY = -0.1, a = 0.1, and eV SD /KT = 0.1. (a) Linear conductance, (b) 
Tunneling current. 



B. Electrochemical potential of voltage probe in weak inter-dot coherent coupling 



Here we discuss the electrochemical potential of the voltage probe ip, which can be measured experimentally. Under 
the mirror symmetry condition, using Eqs. (TT51) . (fTTJj) . (p?5j) . and (J^U) we can show that 



= ^ 2) (0)- 



Furthermore, with the weak inter-dot coupling regime, Eqs. (|18l) and (|19p arc 



$H<t>) 



2HTat c 
4e d 2 + (fiT + hT v ) 2 



Sill ( 



2td 2{16e d 4 - 48e rf 2 ftr y (ftr + HT V ) - (fiT - 3ftr y )(ftr + hT v f}at 

fen i fcT^ \9 ( A - 9. , /t-n i fcT^ \9.13 
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(38) 

(39) 
(40) 



(2) 

From these results, we find that pL v vanishes at cf> = in the linear transport regime and f\i v makes a contribution 

that is independent of the phase </> and inter-dot coherent coupling t c and ahT v . fx^ and /Mp are observables and 
their AB oscillations can be easily detected in experiments. In particular, it is interesting to note that the sign of the 
1st term in Eq. (|40[) depends only on the QD energy e^- Therefore, the sign of the phase-independent contribution of 

(2) 

IMp (4>) is positive (negative) when the QD energy is above (below) the Fermi level. 
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C. General properties of nonlinear conductance coefficient 

Here we discuss more general properties of the lowest-order nonlinear conductance coefficient in Eq. (|14l) . 
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As shown in Appendix [C] under the mirror symmetry condition, we find that the lowest-order nonlinear conductance 
coefficient is asymmetric with respect to the flux 
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Therefore, the lowest-order nonlinear conductance coefficient has no symmetric component. This contrasts with 
the lowest-order nonlinear conductance coefficient in two-terminal system o 18 ' 41 and Mach-Zchnder interferometers^ 2 -, 
which has symmetric and antisymmetric components. 

The position of the current peak or dip at zero magnetic field shifts since the phase symmetry is broken. To 
examine the direction of such a phaseshift, we estimate the lst-order differential coefficient for the lowest-order 
nonlinear conductance coefficient at = 
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In general, we have 
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For clarity, we consider the mirror symmetry and thus obtain 



dG% D {4>) 



d<b 



(fir) 3 (fir, 



= 



h |fi 2 A(e = O,0 = O)| 



-t c a{t c a - e d ). 



(45) 



From Eq. (|4"5"|) . the factor i c a(i c a — e<j) determines the sign of 



Then, we define this factor as 5(a). 



6=0 



In a recent experiment, it was reported that the sign of the coherent indirect coupling parameter a can be changed 
by tuning the gate voltage^!. In the following, as an example, we investigate the direction of the phaseshift when we 
change the sign of a. Then we estimate 



5(a)5(-a) = (t c a) 2 [(i c a) 2 -e d 2 



(46) 



When \t c a\ < \ed\, the slope of the lowest-order nonlinear conductance coefficient at = depends on the sign of a. 
As an example, we show the AB oscillations of the lowest-order nonlinear conductance coefficient in Fig. [5] (a) when 
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FIG. 5: AB oscillations of the lowest-order nonlinear conductance coefficient when Ts = T_d = = T, t c /hT 
ea/hT — 1. The solid and broken lines indicate cases where a = 0.5 and a — —0.5, respectively. 



-1, and 



0.02 




0.18 



8, 

-5 




(a) 



(b) 



FIG. 6: AB oscillations of the lowest-order nonlinear conductance coefficient when Ts = T_d = = V, t c /hT = — 1, and 
€d/KT = 0.2. The solid and broken lines indicate cases where a = 0.5 and ot — —0.5, respectively. 



Ts = I\d = r v , t c /hT = —1, \a\ = 0.5, and €d/KT = 1. In contrast, when |t c a| > | | , the slope of the lowest-order 
nonlinear conductance coefficient at <f> = is independent of the sign of a. As an example, we plot the AB oscillations 
of the lowest-order nonlinear conductance coefficient in Fig. [6] (a) when Ts = Tp = T v = T, t c /HT = —1, \a\ = 0.5, 
and e d /hT = 0.2. 

From the above results, we consider the direction of the phaseshift for the AB oscillations in the current through 
a DQD. In Fig. [5] (b), we plot the AB oscillations in the current through a DQD when ed/hT = 1, t c /hT = — 1, 
Ts = Td = T v = r, and \a\ = 0.5, namely \t c a\ < \ed\- For a = 0.5, the conductance exhibits a dip at <p — as 
shown in Fig. [5] (a). According to Fig. [5J the position of this dip shifts to the negative phase direction in a nonlinear 
transport regime. Similarly, for a = —0.5, we have the conductance peak at <p = as shown in Fig. [2] (a). According 
to Fig. 0(a), the position of this peak should shift to the negative phase direction in a nonlinear transport regime. 
The direction of the phaseshift for the current through a DQD at eVsn/^T = 2 is consistent with a prediction for the 
lowest-order nonlinear conductance coefficient as shown in Fig. [5](b). 

Next we consider the situation when ed/KT = 0.2, t c /HT = —1, Ts = Tjj = T v = T, and \a\ = 0.5, namely 
\t c a\ > \ed\- For both a = 0.5 and a = —0.5, the conductance shows a dip at 4> = as shown in Fig. [2](b). According 
to Fig. [6] (a) , the position of these dips should shift to the negative phase direction in the nonlinear transport regime. 
The direction of the phaseshift for the current through a DQD at eVso/KT = 2 is consistent with a prediction for the 
lowest-order nonlinear conductance coefficient as shown in Fig. [6] (a). 



From Eq. (|45|) . when t c a = t d , we have 



= 0. However, the lowest-order nonlinear conductance 

4>=o 



11 



-e- 

si 




I -0.01 



FIG. 7: Behavior of the lowest-order nonlinear conductance coefficient under the mirror symmetry condition when t c a = 
(td/fiT = —0.5, t c /hY = — 1, a = 0.5, and T v = F). </> = is an inflection point, and the lowest-order nonlinear conductance 
coefficient has an antisymmetricity under the mirror symmetry condition (ei = £2 = €d, Ts = Yd = Y, and Y^i = Y V 2 = Y<p). 



coefficient does not have an extreme value at <p = since we have 



= under the same condition. As 



= 



a result, = is an inflection point as shown in Fig. [7] when td/KY = —0.5, t c /HT = — 1, a = 0.5, and r v = T. 



V. DEPHASING EFFECTS OF VOLTAGE PROBE ip 



So far we have discussed the phaseshift in the nonlinear transport regime. In this section, we consider the amplitude 
of the AB oscillations in the transport properties to study dephasing effects induced by the voltage probe. To observe 
the AB oscillation, we need the coupling T v . However, this causes dephasing of the electronic states in the DQD. In 
this section, we study the competition between these two antithetic concepts. In particular, we examine the a and T v 
dependences of the AB oscillations in the linear and nonlinear conductances. For simplicity, we focus on the mirror 
symmetry in this section. 



A. Linear conductance 



Here we discuss the AB oscillations in the linear conductance. To investigate the interplay between the two antithetic 
concepts as mentioned above, we examine how the coherence is modulated as the coupling T v increases. As a physical 
quantity that characterizes the coherence, we define the visibility of the AB oscillation as follows 

„ _ u DQD,max u DQD ,min , 
DQD,max ' u DQD,min 

where G^q D max and G^q D min correspond to the maximum and minimum values in the AB oscillations of the linear 
conductance, respectively. In Fig. [8j we show the coupling r v dependences of the visibility of the AB oscillations in 
the linear conductance for various \a\ values when ed/fiT = —t c /hT = 1. In the weak coupling regime, the visibility 
increases monotonically as the coupling T v increases. This result reveals that we need a stronger coupling T v to 
observe the AB oscillations with higher visibility. However, for \a\ ^ 1, the visibility has a maximum value, and 
the visibility decreases with T v in the strong coupling regime. This result means that the strong coupling gives 
rise to the dephasing, and leads to the loss of the coherence. The interplay between these two features provides the 
maximum visibility as shown in Fig. [U Moreover, in the limit of T v ^> T, the visibility of the AB oscillations has the 
following leading term in an asymptotic series for \a\ 1 

8\e d t c a\{3 + a 2 ) ( 1 ' 



{I -a 2 ) 2 \hT 



Vr^oo- ' " ' ttH • (48) 



Vr^-yoo decreases with (r v ) 2 with a monotonically increasing coefficient \a\ 
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FIG. 8: Visibility of the AB oscillations in the linear conductance as a function of Y v for various \a\ values when td/KT 
-to/hT = 1. 



In contrast, for |cv| = 1, there is no maximum visibility, and the visibility increases monotonically as the coupliing 
T v increases. In the limit of infinite T v , the visibility of the AB oscillations with |cv| = 1 can be expressed as 



h 






i = 7T)) 


e 2 
h 


+ min {G^qjj^ 




i = 7T)) 



where 



^-°= = hTT^- — 777" "V- ( 49 ) 



^-0)=' T _^ (50) 



G DQD (4> = *) - ^ 4{ed + tc y + {hrr (51) 



In the situation shown in Fig. we have Vr -+oo = 8/9. Moreover, when T ^ 0, using the relation 

^ ifawbffF - 1 ' (52) 



we can prove that the visibility Vr^-^oo < 1 from Eq. The derivation of the expression of the visibility in Eq. 

is given in Appendix [Dj Such exceptional behavior at \a\ = 1 can be explained as follows. With the tunnel-coupled 
symmetric and antisymmetric states as a basis, the linewidth function matrices are given by 

r s = U] ]) , r D = U 1 -i) , T*{<t>) = rJ 1 + acos t / asin M. (53) 

2 y 1 1 y 2 \ — 1 1 J \ — lasmcp 1 - a cos <p J v ' 

Then, the coupling strength between the symmetric (antisymmetric) state and the voltage probe is characterized by 
r v (l + a cos 4>) (r^(l — a cos </>)). Thus, at <f> = (<f> = 7r), the antisymmetric (symmetric) state is dephasing-free from 
the voltage probe, and the coupling T v helps to enhance the coherence. As a result, when + t c = (e<j — i c = 0), 
the resonant tunneling process is realized through the symmetric (antisymmetric) state in the series-coupled DQD, 
and the linear conductance has a value of e 2 /h. This is the origin of the high visibility in the limit of \a\ = 1 and 

— > co. In contrast, for \a\ ^ 1, both the symmetric and antisymmetric states are dephased by the coupling 
with the voltage probe. Consequently, the linear conductance vanishes without depending on 0, and the visibility of 
the AB oscillations becomes zero. 



B. Lowest-order nonlinear conductance coefficient 

(2) 

Here we discuss the lowest-order nonlinear conductance coefficient G D q D {4>) with respect to the source-drain bias 
voltage. In Fig. H](a), we plot the AB oscillations of G D Q D ((j>) for ed/hT = 1, t c /HT = —1 and a = 0.5 for various 
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FIG. 9: AB oscillation of the lowest-order nonlinear conductance coefficient G^q D (0) and T v and the \a\ dependences of its 



amplitude when e^/KT — 1, t c /hT = —1. (a) AB oscillations of G^Q D {(f)) for a = 0.5. (b) T v dependences of AG^q D for 



-.(2) 



various \a\ values, (c) \a\ dependence of the maximal value AG^q D max in (b) 



The amplitude of the AB oscillation depends on T v values. We discuss the T^, dependences for the amplitude of the 
AB oscillations in G D q D (4>), which is defined as 

^Gdqd = GpQ D (</) max ), (54) 

(2) (2) 

where <ft ma x is the phase when G D q D (4>) is maximal in the AB oscillation. When increases, AG D q D has a peak as 
shown in Fig. E] (b) when e^/fiF = 1 and t c /hT = — 1. To understand the behavior of G^q£,(0) in the weak coupling 
regime (T v <C T) in Fig. |H] (b), we consider the situation where e,i = —t c , and \a\ <C 1. Then, we obtain 



32i c 2 |sin( 



128i c 4 + 2At c 2 (hT) 2 + (hT) A 



Up to the order of |a|, we have 



_/ § (a>0) 
-f (a<0) 



(55) 



(56) 



(2) 

Thus, the slope of AG D q D for T v <C T increases with a|. 

In the limit of r v 3> T, the amplitude of AB oscillations in the lowest-order nonlinear conductance coefficient has 
the following leading term of an asymptotic series for \a\ ^ 1 

A ^,(2) e 3 32(fiF) 2 \t c a [e d - i c cv cos(0 max )] sin(^ max )| / 1 \ 5 

AG ^ ~ x o^ 5 ? { ) ■ (57) 

When = — t c and |a| <C 1, we have 

(2) e 3 32(fiT)V|q| 
AG DQD ~ - (fir ^ )5 . (58) 

Under the condition where |a| C 1 is not satisfied, it is difficult to discuss the asymptotic behavior of AG^q D 

(2) (2) 

since </> ma x is a function of a and T^. In general, AG D q D is a function of \a\ since G D q D is invariant under the 
transformation that we change the sign of a and shift the phase by ir. Moreover, we plot the |a| dependence of the 
peak height of AG^q D (indicated as AG^q D max ) as shown in Fig. HJc). The peak height increases monotonically 
as | a | increases in the same way the visibility of the AB oscillation in G^q D {4>). 

VI. SUMMARY 



We studied linear and nonlinear three-terminal transport through an AB interferometer containing a DQD using the 
nonequilibrium Green's function method. We introduced coherent indirect coupling between two quantum dots via a 
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voltage probe ip. The linear conductance exhibits phase symmetry without depending on the various parameters of the 
model. However, in the nonlinear transport regime, the phase symmetry is broken and the phase of the AB oscillations 
shifts. We showed that the lowest-order nonlinear conductance coefficient with respect to the bias voltage contributes 
to the phascshift. In particular, when jal^r^, \t c \ <C hT, where we can neglect the higher harmonic components of the 
AB oscillations, we proved that the sign of the lowest-order nonlinear conductance coefficient is directly related to that 
of the phaseshift, and the value of the phaseshift is determined by the quotient between the linear conductance and the 
lowest-order nonlinear conductance coefficient. In the weak coherent indirect coupling and low bias voltage regimes, 
the phascshift is independent of the coherent indirect coupling parameter and monotonically increasing function with 
respect to the source-drain bias voltage. Moreover, we obtained a condition where the direction the phase of the 
AB oscillation shifts from the lowest-order nonlinear conductance coefficient. In the coupling r v dependence of the 
visibility of the AB oscillations in the linear conductance, we found that the visibility has a maximum value except 
when a = l. When a = 1, the visibility is a monotonically increasing function of T^. 
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Appendix A: Derivation of probability conservation 

In this appendix, we derive the relation of the probability conservation given by Eq. (|27|) . The left-hand-side of 
Eq. ([271) is 

$>„ e (e,0) = ^Tr{G r (e,0)r«(0)G a (e,^(0)} 

= Tr {G r (e, </>)r(0)G Q (e, 0)1^(0)} - Tr {G r (e, </>)r"(0)G Q (e, </>)r^)} . (Al) 
Here the Dyson's equation for the retarded Green's function is 

G r M) = g r (e)+g r (e)X r ^)G r (e,ct>) 

= g r (e)-±g r {t)r(4>)Gr(e,<j>), (A2) 

where g r (e) is the retarded Green's function of an isolated DQD. Thus, we have 

r(0) = *{£''(<« -£*(</>)} 

= i {[G a (e,0)]- 1 -[G-(e,0)]- 1 }. (A3) 

Therefore, Eq. (fMJ) is 

5>* C M) = iTr{Cr{e,4) [[G a (e, cf,)}- 1 - [G r (e, cb)}- 1 ] G a (e, 4>)T" (4>)} — Tr {G r (e, <p)T v (<p)G a (e, <fi)T u (<fi)} 

= iTr {G r (e, 0)1^(0) - G a (e, <f)T" {<)>)} - Tr {G r (e, 0)r^(0)G a (e, <t>)T" (</>)} 
= *Tr {r^)GXe, <f>) - r^)G a (e, <f>)} - Tr {G r (e, W(<f>)G a (e 7 W (</>)} 

= zTr {r"(0)G°(e, <t>) ([G a (e, cf)]- 1 - [G r (e, 0)]- 1 ) G r (e, </>)} - Tr {G r (e, <f)T» {4>)G a {e, W(<f>)} 
= Tr{r-(0)G a ( e ,^)r(0)G-( e ,0)} -TriG^e^r^G^^r^)} 
= ^TrjG'^^^G^e,^^)} 

= ^T ?I/ (e,# (A4) 
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Appendix B: Derivation of phase symmetry 



In this appendix, we prove the phase symmetry relation ([28]) . Using Eqs. ([20]) and ([25]) , the linear conduct 



. e 2 r SD (e = 0, 0)T yS ( £ = 0,0)+ T S p(e = 0, 0)?>(e = 0, 0) + T Sy (e = 0, 0)T yD ( 



e = 0, 



1^(6 = 0,0) + 1^ = 

Using relation (f2"T|) , the numerator of the right-hand-side of Eq. (|B1I) can be rewritten as 

T S d(£ = 0, 0)T vS (e = 0, 0) + [T SD (e = 0, 0) + T Sv {e = 0, 0)] T^(e = 
= T SD (e = 0, <j))T vS {e = 0, 0) + [T DS {e = 0, 0) + T vS {e = 0, 0)] T vD (e = 

By continuous use of Eq. ([27]), Eq. (|B2]l is 

[r S u(e = 0, 0) + l>(e = 0, 0)] T^e = 0, 0) + T DS (e - 0, 0)-7>(e = 
= PWe = 0, 0) + T Dv (e = 0, 0)] T vS (e = 0, 0) + T DS (e - 0, 0)2>(e = 
= T DS (e = 0, 0) [T vS (e = 0, 0) + 7>(e = 0, 0)] + T^e = 0, 0)T^(e = 

Therefore, by comparison with Eq. (|20[) . we obtain the phase symmetry relation 

G% D {-<t>) = G% D {4>). 



;ance is 



(Bl) 



(B2) 



(B3) 
(B4) 



Appendix C: Derivation of antisymmetricity for lowest-order nonlinear conductance coefficient 

Here wc show that the lowest-order nonlinear conductance coefficient is asymmetric with respect to the flux. Using 
relations (|25|) and (|26j) . when the direction of the magnetic flux is reversed, the flux-dependent contribution in the 
1st line of the right-hand side in Eq. (|4"T]) is 



Tsv(e = 0, 



T Dip (e = 0, 



T Dv (e = 0, 0) - T Sip {e = 0, . 



(CI) 



T Slp {e = Q,-<t>)+T Dlp (e = Q,-<i>) J I T Dlp (e = 0, 0) + T Slf {e = 0, 
Then, this contribution is symmetric with respect to 0. Similarly, the 2nd line of the right-hand side in Eq. (|41[) is 



dT vS {e, 



ck 



dT Sv (e,-<P) 




dTn v (t,-(f>) 




de 


e=0 


de 


e=0 



=0 Ts<p(e = 0,-4>)+T Dv> (e = 0, 



T, 



0,-0) 



de 



dT vS (e,<t>) 








de 


e=0 


de 


e=0 





e=0 



=0 7Vs(e = 0,0) +T Sv (e 



0. 



Ts v (e = O,0) 



de 



e=0 



T (pS (e = 0,(f>)+Ts v (e = 0,4>) 

[T< ?p (e = O,0)+T /3¥ ,(e = O,0)] + 



dT Dv ,{e, 



de 



T v s(e = 0, 



c9r 



T^(e = O,0)+T Sv ,(e = O,0) 
■?Vs(e = 0,. 





^ 


dTn v (e,4>) 




de 


£=0 


de 


e=0 



e=0 



T Sv (e = O,0)+T D¥ ,(e = O, l 



(C2) 



Then, this contribution is asymmetric with respect to 0. As a result, from Eqs. (jUJ), (|C1I) and (|C2[) . we find that 



-G 



(2) 
DQD 



Appendix D: Derivation of visibility (|49[) 



In this appendix, we derive expression ([49]) for the visibility of the AB oscillations in the linear conductance in the 
limit of a — 1 and T v — > oo. We consider the mirror symmetry (ei = £2 = e<j, r<j = = T, and T^i = T V 2 = Ty)- 



Under this condition, the linear conductance is given by 



.,2 



G'd'qdW = -^T DS (e = 0,, 



(hTf 



h 4(e d -t c cos(/)) 2 + (^r) 2 



since, in the limit of a = 1 and — > oo, we have 



and thus the second term in Eq. ([20]) vanishes. Then, from the following condition 



we have the condition of d> for extreme values: 



Therefore, from the definition (1471) , we obtain the expression of visibility (|49 
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(Dl) 



T vS (e = O,0) = T Dv (e = O,^»)<x-^ ! (D2) 

i v 



— of = °> ( D3 ) 



sin cj) = , cos</> = — . (D4) 
If e c i ^ \t c \, the condition cos^) = Cd/t c leads to </> ^ nir, where n is an integer. Then, we have 

2 

G% D . max = \ , G% D>min = min (G% D {4> = 0), G% D {cj> = tt)) . (D5) 
Similarly, for = \t c \ 1 the condition cos</> = £d/t c leads to <f> = or ir. Then, wc have 

G% D , max = £ , G«^ mm = min (g% d {4> = 0), G% D (<f> = tt)) . (D6) 
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